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Motivated by the reported peculiar dynamics of a red blood cell in shear flow, we develop 
an analytical theory for the motion of a nearly-spherical fluid particle enclosed by a 
visco-elastic incompressible interface in linear flows. The analysis explains the effect of 
particle deformability on the transition from tumbling to swinging as the shear rate 
increases. Near the transition, intermittent behavior is predicted only if the particle has 
a fixed shape; the intermittency disappears for a deformable particle. Comparison with 
available phenomenological models based on the fixed shape assumption highlights their 
physical foundations and limitations. 
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1. Introduction 

A membrane-enclosed particle, such as the red blood cell (RBC), exhibits rich dynam- 



ics in flow (see for recent reviews Abkarian & Viallat (2008); Vlahovska et al. (20096); 



|Guido fc Tomaiuolo (20091). Two classic types of behavior in steady shear flow are 



(1) tank-treading (TT), in which the RBC shape is steady and the membrane rotates as 
a tank-tread; the cell major axis is tilted with respect to the flow direction and the incli- 
nation angle remains fixed in time, and (2) tumbling (TB), in which the RBC undergoes 
a periodic flipping motion. Recently, a new type of motion called swinging (SW) has been 
experimentally observed ( Abkarian et aZ.|2 007). In this case, the RBC's tank-treading is 
accompanied by small oscillations in the inclination angle. Similar phenomenon has been 



reported for capsules (Walter et al. 2001[ Ramanjuan fc Pozrikidis| [l998| |Navot| |1998l 
and drops covered with adsorbed protein layer ( Erni et q/.|2005 ) 



The variety of RBC's motions stems from the unique mechanical properties of the 
interface. The cell membrane is made of a lipid bilayer attached to an underlying spectrin 
network. The lipid bilayer behaves as a two-dimensional incompressible fluid, while the 
polymer network endows the membrane with shear elasticity. Closed lipid membranes 



(vesicles) display TT, TB but no SW (Kantsler fc Steinberg 2005 Mader et al. 2006 



Kantsler fc Steinberg|2006 Deschamps et a/.|2009a|6|). The inclination angle can oscillate 



around the flow direction (Kantsler & Steinberg 2006 Mader et aL|2006 Barthes-Biesel 
& Sgaier 1985), but this breathing (VB, also called trembling) motion differs from the 
swinging observed with RBCs; in the latter case, swinging does not necessarily involve 
shape deformation. 

The type of motion a vesicle or RBC undergoes depends on the viscosity mismatch 
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between the inner and suspending fluids, applied shear rate, and the non-sphericity of 
the rest shape. For a given viscosity ratio and shape, at low shear rates the RBC tumbles 
because elastic tensions immobilize the interface causing the cell to behave as a solid 
object. At high shear rates the applied stress overcomes the elastic tensions and drags 
the membrane in motion. As a result, the RBC adopts a steady TT shape whose major 
axis "swings" . This motion originates from the unstressed non-spherical shape of the 
RBC, in which membrane elements at the equator and the poles are not equivalent. 
During one TT-period, an element passes twice through its unstressed position where it 
releases elastic energy. As the shear rate increases, the amplitude of the angle oscillations 
decreases. 

Phenomenological studies of swinging (Abkarian et q/.|2007 Skotheim & Secomb 2007 



Kessler et aL|2009| |Noguchi|2009| based on the theory by |Keiier fc Skalak| ( |1982| ), which 



models the RBC as an ellipsoid of fixed shape, qualitatively capture the physics of the 
phenomenon. However, the quantitative predictions match poorly with simulations and 



experiments. For example, the model proposed by Abkarian et al. (2007) agreed with 



the experimental data only if the shear elastic modulus of the membrane was assumed 



to be lower than the generally accepted value. Furthermore, Skotheim & Secomb (20071 



predicted intermittent behavior, for which no evidence was found in the numerical sim- 
ulations ( ]Kessler et oT1|2008| |Sui et q/.||2008||Bagchi fc Kalluri||2009[ ) . A r ecent detailed 
analysis constructed the dynamical phase diagram of the reduced model ( Kessler et al. 



2009) in both steady and oscillatory shear. 



The purpose of this work is develop a rigorous analytical theory that disposes of the as- 
sumptions for fixed ellipsoidal shape and compressible membrane inherent to the studies 
based on the Keller-Skalak model. We generalize the theory for the dynamics of a vesicle 
made of fluid incompressible membranes ( Misbah||2006| Vlahovska fc Gracia||2007] ) to 
include membrane shear elasticity. Section |2| formulates the model, Section |3j summarizes 
the theory for fluid membranes (vesicles) , Section El discusses the evolution equations 
for the particle shape and orientation angle, and Section [5] analyzes the effects of shear 
elasticity on particle dynamics. 



2. Problem formulation 

We model the RBC as a closed membrane ( "capsule" ) with total area A. The membrane 
encapsulates a fluid of viscosity Xrj and it is suspended in a fluid of viscosity 77; A denotes 
the viscosity ratio. Both interior and exterior fluids are incompressible and Newtonian. 
The particle has a characteristic size Rq defined by the radius of a sphere of the same 
volume. The nonsphericity of the rest shape is characterized by an excess area 



A = A/Rl - 4tt . 

The typical value for a RBC is A ~ 4. 

The RBC is placed in a steady two-dimensional linear flow 

v°°(r) =7J/-x, 

where 7 is the strain rate. A sketch of the problem is shown in Figure [T] 



(2.1) 



(2.2) 



2.1. Fluid motion and fluid-membrane coupling 

At the length scale of the micron-size RBC, water is effectively very viscous and creeping- 
flow conditions prevail. Fluid velocity v( Q ) and pressure p^ of the interior ("a = in") 
and suspending ("a = ex") fluids obey the Stokes equations and the incompressibility 
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v = 72/x 

Figure 1 . A. Sketch of the streamlines of a linear shear flow. The shear flow is a superposition 
of pure straining flow and rigid body rotation B. A soft particle deforms to an ellipsoid in a 
simple shear flow that can tank-tread (ip = const) or tumble (ip continuously increases) 



condition 

V-T (q) =0, V-v (q) =0, 
where T is the bulk hydrodynamic stress 
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(2.4) 



I denotes the unit tensor and the superscript f denotes transpose. Far away from the 
particle, the flow field tends to the unperturbed external flow v cx —> v°°. The velocity 
field is continuous across the membrane. However, the hydrodynamic stresses undergo a 
jump, which is balanced by membrane surface forces 



n 



= t B 



(2.5) 



where n is the outward unit normal vector and the membrane surface forces t mm arc 
discussed next. 

2.2. Mechanics of flexible fluid membranes 

The RBC membrane consists of a lipid bilayer attached to protein scaffold. The lipid 
bilayer endows the membrane with fluidity, incompressibility and bending rigidity, while 
the polymer network gives rise to resistance to shearing. The membrane thickness is 
5nm, or about 1/1000 of the cell radius. Due to the large separation of lengthscales, the 
membrane can be treated as a two-dimensional surface embedded in a three-dimensional 
space. 



Within the framework of the minimal model (Seifert 1997), the bending resistance 
gives rise to a surface force density 



-K (AH 3 - AKH + 2V 2 S H) n , 



(2.6) 



where H = (1/2) V • n and K = (1/2) [(V • n) + Vn : Vn^] are the mean and Gaussian 
curvatures, and k is the bending modulus. The surface gradient operator is defined as 
V s = I s ■ V, where the matrix I s = I — nn represents a surface projection. 

The bilayer consist of fixed number of lipids, which are optimally packed with fixed 
area per lipid (under moderate stresses). As a result, an element of the bilayer membrane 
only deforms but can not change its area. Under stress, the membrane develops tension (a 
two-dimensional pressure) , which adapts itself to the forces exerted on the membrane in 
order to keep the local and total area constant (Seifert 1999). This tension is non- uniform 
along the interface and varies with forcing. The corresponding surface force density is 



2aHn- V s a. 



(2.7) 



where a denotes the local membrane tension. 

For lipid bilayers in the fluid phase, the lipids are free to move within the monolayer. 
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Therefore, a pure fluid bilayer membrane is infinitely shearable (Dimova et al 



2006). 



However, the polymer network lining the bilayer develops elastic tensions when sheared. 
Various constitutive laws exist to model the elastic behavior of RBCs and capsule mem- 
branes (Pozrikidis 2003 Barthes-Biesel||1991[). Assumin g a linear elastic behavior, the 
elastic tractions are given by ( Barthes-Biesel fc Rallison|1981| [Edwards et aL||1991[) 



t" = - 2(K A - M )(V S • d)Hn + (K A - ^)V S V S • d + M V S ■ [v s d • I s + I s • (V s d) f 

(2.8) 

where d is the displacement of a material particle of the membrane. Ka is the stretch 
and n is the shear elastic moduli. For RBC n ~ 10 -6 N/m and Ka ~ 2Q0N/m (Dimova 



et al. 2006 
& Rallison 



We can identify u-i — Ka — /U > a 3 — i 1 m the notation of |Barthes-Biesel 
(11981). In general, a in Eq. (2.7) includes both tensions arising from shear 



elasticity and area-incompressibility. For an area-incompressible membrane the first two 
terms in Eq. (2.8) vanish because V s • d = 0. 



2.3. Time scales 

Viscous forces exerted by the extensional component of the flow act to distort the shape 
on a time scale 

7} = (1 + X)r l • (2.9) 

Several intrinsic relaxation mechanisms oppose the deformation. Bending stresses work 
to bring the shape back to its preferred curvature state; the corresponding time scale is 

W* . (2.10) 



1 



(2.11) 



Relaxation driven by shear elasticity occurs on a time scale 

(1 + A)?/i? 
T v> = • 

The factor (1 + A) reflects the fact that the more viscous fluid controls the dynamics. In 
shear flow, particle rotation away from the extensional axis of the imposed flow effectively 
decreases the effect of the straining; the associated time scale is 



7 



(2.12) 



The strength of the relaxation mechanisms that limit shape deformation by the flow is 
quantified by the corresponding dimensionless parameters: the capillary number 



Ca Kl = 



(2.13) 



elastic capillary number 



and the rotation parameter 






£ = V 



1 



(2.14) 



(2.15) 



The interplay of these time scales leads to the complex dynamics of RBCs and vesicles. 
The elastic capillary number x _1 plays the role of a dimensionless shear rate. For a fluid 
vesicle x = an d the dynamics is relatively insensitive to the shear rate ( Misbah||2006 ); 
for a given excess area, the TT-TB transition of a vesicle is controlled by the viscosity 
ratio A. 



Swinging of red blood cells in simple shear 



Henceforth, all quantities are non-dimensionalized using n, Rq, and 7. Accordingly, the 
time scale is 7 -1 , the velocity scale is jRo, bulk stresses are scaled with 777. 

2.4. Perturbative solution for a nearly-spherical shape 
In order to solve the problem analytically, we consider a nearly-spherical particle shape, 
i.e., A <C 1. In a coordinate system centered at the particle, the radial position r s of the 
interface can be represented as 



r B = l + /(M,*) ; 



(2.16) 



where / is the deviation of particle shape from a sphere. The exact position of the interface 
is replaced by the surface of a sphere of equivalent volume, r — 1, and all quantities that 
are to be evaluated at the interface of the deformed particle are approximated using a 
Taylor series expansion. The leading order analysis for fluid vesicles has been done in 



Vlahovska & Gracia (2007). Here, we modify the solution to include shear elasticity. 



In Eq. (2.161, the function / representing the shape depends only on angular coordi- 



nates. Thus, it is expanded into series of scalar spherical harmonics Yj m (El) 



J — 2^i / , JjmYjm , (2.1 7) 

In the above equation, the summation starts from nonzero j '■ = 2 because j = 1 correspond 
to translation of the center of mass. The constraint on fixed total area serves to relate 
the amplitude of the perturbation / and the excess area A 

(2.18) 



A = / ^dn -4n = J2 °' + 2) 2 ' 1] f 3m fU + 0(f) 

jm 



where f denotes the unit radial vector and the sum over j starts from 2, \m\ < j and 
/*,„ = (— l) m fj- m . Since the rest shape of the particle is characterized by small excess 



1 jm 

area 



/ ~ A 1 / 2 « 1. 



Dynamics of a nearly spherical closed fluid membrane in a shear 
flow 



The theory for vesicles, \ = 0j is well developed (Misbah 2006 



2007 Lebedev et al. 2008 Danker et al. 2007 Kaoui et al. 2009 



Vlahovska & Gracia 



Schwalbe et al. 2010 ) 



Here we summarize the main results. In shear flow, a vesicle deforms into an ellipsoid and, 
hence, its shape is specified by the j ' = 2 spherical harmonics. In the flow plane x — y, the 
vesicle shape / is characterized by two components, /2±2 corresponding to deformation 
along the flow axis x and the straining axis x — y. The out-of-plane deformation along the 
vorticity z axis is described by the /20 mode. For simplicity, fo±i modes are neglected. 
Instead of shape modes, the vesicle dynamics can be more conveniently described in 
terms of the orientation angle, tp, and R, which measures the ellipticity of the vesicle 
contour in the x — y plane ( |Misbah||2006[ ) 

/ 2±2 = i?exp(T2k/>). (3.1) 



The /20 is slaved to the / 2 ±2 modes and it is determined from the area constraint (J2.18 ) 



/20 — 



A 



— 2/22/2-2 



1/2 


FA 


- 


= 




-2R 2 




['2 





1/2 



(3.2) 
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The evolution equations for the shape and orientation of a fluid membrane vesicle in a 
simple shear flow are (Misbah 2006) 

dip _ 1 h 



2R(t) 



cos [2ip{t)} 



(3.3) 



OR 
~dt 



= h 



I - 4^- ) sin pty(t)] 



(3.4) 



Unlike drops and initially spherical capsules, the vesicle 



where h = 

motion described by equations (|3.3| and (3.4 1 is nonlinear (even at leading order) and in- 



dependent of the elastic properties of the interface. This is due to fact that the interfacial 
stresses are dominated by the tension that arises from the surface-incompressibility. 

This set of coupled nonlinear equations has a stable fixed point corresponding to the 
tank-treading state (R* = \fK/2 , ip* = arccos y/A/h) and a closed orbit centered at 
(ip* = , R * = h) describing the breathing mode. Tumbling corresponds to no equilibrium 
point. The TT fixed point loses stability at a critical viscosity ratio 



A, = -: 



32 
23 



120 
~23 



2?r 
15A 



(3.5) 



If there is no deformation along the vorticity direction, i.e., /20 = at all times, Eq. 
(3.2) implies that R remains constant and equal to its maximum value vA/2. This 



situation resembles the Keller-Skalak model: the vesicle shape is a fixed ellipsoid and the 
vesicle dynamics is described only by the variations of the angle ip (note, however, that 
unlike the Keller-Skalak solution, our velocity field is strictly area-incompressible). In 
this case the bifurcation is from TT to TB. 

If /20 7^ 0: the transition is from TT to VB. In the breathing (VB) mode, the vesicle 
undergoes periodic shape deformations along the vorticity direction. As a result, the 
vesicle appears to tremble in the flow direction. 



As already discussed by Lebedev et al. (2008) and Kaoui et al. (2009), within the 



leading order theory ( 3.3 ) and ( 3.4 ) the TB and VB modes coexist, and the mode selection 



is determined by the initial conditions. 



4. Effect of the shear elasticity 

The inclusion of solid-elastic deformation in the fluido-dynamic problem is, in general, 



not a trivial task (Barthes-Biesel 1980). Since fluids have no memory, their motion is 



described in a fixed laboratory frame (Eulerian approach). However, membrane elastic 
stresses depend on the displacement of material particles, i.e., the membrane deformation 
has to be, in principle, described in a material (Lagrangian) frame. For example, at steady 
tank-treading state, the deformation at a fixed Eulerian point is constant in time unlike 
that of a fixed material point because material elements rotate. Coupling the membrane 
deformation and fluid motion requires transformation between the Lagrangian and the 
Eulerian representations. A nice discussion of this problem is given by Barthes-Bie sel fc| 
Rallison (1981); we closely follow their approach. 



4.1. Reference unstressed shape 

The position of a membrane element in the unstressed configuration is labeled by X. At 
time t, its position in the Eulerian frame is x(X,£). The element displacement is 

d = x(r,M,t)-X(r,M,t) (4-1) 
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The motion of the interface is parametrized with respect to a sphere 

d(t) = d t + r s {t)r . (4.2) 

where d t is tangential to a sphere. Since the interface is area-incompressible, the radial 
and tangential displacement components are not independent. Accordingly, the radial 
displacement is sufficient to describe the motion of a material point; the tangential dis- 
placement is determined from the condition V s ■ d = 0. The unstressed shape is assumed 
to coincide with the initial shape and it is specified by 

r ref {t) = 1 + g jm Y jm (e, cf>) . (4.3) 

The imposed shear flow consists of straining flow and rigid body rotation, which does 
not generate deformation. Accordingly, the distortion of a material element needs to be 
defined relative to a rotating unstressed configuration. At leading order, the reference 
unstressed configuration rotates with the flow vorticity 

dg 2m .m 

~dT = 1 ¥ 52m ' (4 ' 4) 

4.2. Evolution eguations 



The solution follows standard steps (Vlahovska fc Gracia||2007) and it is summarized in 



Appendix [Dl We obtain that the shear elasticity modifies the evolution equations (3.3 1 



and (3.4 1 as follows 



S = 1 r- cos \2ib(t)] + x T -7= sin \2<b(t) - 2ib(t)} (4 5) 

dt 2 2R(t) LVWJ X R(t)V30n L ' WJ l ' 



9R_,f, Mt)' 



;jj -h^l~4^ym[2yj(t)}+ X ^^{r (A-AR(t) 2 )cos[2ct>(t)-2^(t)} 

-2R(t)f 20 (t)g%} 

where 



f 20 (t) = VA/2 - 2R{tf (4.7) 

and ro 7^ defines a non-spherical reference shape 



g° 2±2 = r exp(T2i^o) .g 2 °o = \/ A / 2 " 2r o ' ( 4 - 8 ) 

where g^ = gj m (t = 0). The tank-treading frequency is 



(4.9) 



9^ _ _1 
~dt ~~2 : 

where </> is the angle between the position vector of a material particle and the flow 
direction. If the reference (unstressed) configuration of the vesicle is a sphere, both r = 
and (?2o = 0. In this case the evolution equations become independent of the shear 
elasticity, similarly to the way bending elasticity scales out from the evolution equations 
for the vesicle dynamics (Olla 2000). 



5. Results and discussion 

In this section we first discuss the reduced model corresponding to a non-deformable 
ellipsoid in order to compare our work to earlier studies based on the Keller-Skalak's 
theory. After that we proceed to analyze the more general case of a deformable particle. 
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(a) 




(b) 




Figure 2. Swinging motion: (a) the inclination angle oscillates in time, (b) the limit cycle in 
the phase portrait. The amplitude of the oscillations increases with the elasticity number x — 1 
(red) , x — 10 (blue) x = 29 (black). The other parameters are: A = 0.02 and A = 5 . 



5.1. Dynamics of a particle with a fixed ellipsoidal shape 



The evolution equations (4.5 1 and (4.6 1 reduce to a shape-preserving model, if there is no 
out-of-the-shear-plane deformations, / 2 o = const. Because the evolution of the /20 mode 



is slaved to the 2 ± 2 modes (as seen from (4.7)), this condition can be strictly enforced 
only if /20 = and 520 = 0- This implies that the ellipticity of the vesicle contour is 
constant 

R(t) = r = VA/2 . (5.1) 

The particle dynamics is described by only one variable, the inclination angle ip 

d^ 1 h 2h 

~dt 



- H — = cos (2t/>) - x- 



.8in(t + 2ip-2ip ) 



(5.2) 



Examples of the angle evolution are shown in Figures 2](6 At high shear rates (small 
x), the particle is swinging; the inclination angle oscillates around a positive value as 
shown in Figure [2Ja. Eq.(5.2) has a stable limit cycle as seen in Figure [2]b. As the shear 
rate decreases (% increases) the amplitude of the oscillations increase and the radius of 
the limit cycle also increases. Eventually the limit cycle loses stability and the particle 
begins to tumble. The tumbling dynamics is illustrated in Figures[3ja and[3]b. Increasing 
the viscosity ratio, at a fixed elasticity, also causes the limit cycle to disappear in favor 
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(a) 



(b) 




%|«"l-5 



Figure 3. Tumbling motion: (a) the inclination angle decreases continuously with time, (b) the 
phase portrait shows the angular velocity as a function of the inclination angle. Close to the 
transition the angular velocity is nonuniform. As the elasticity increases the limiting behavior of 
rigid ellipsoid is approached. Parameters: \ — 30 (red) , % = 100 (blue), and \ — 1000 (black). 
The excess area is A = 0.02 and viscosity ratio is A = 5. 



of tumbling because increasing the viscosity ratio decreases the mean inclination angle. 
However the swinging amplitude is relatively insensitive to the viscosity ratio, as seen in 
Figures [4] a and|4]b. 

At intermediate values of the elasticity parameter and high viscosity ratios, the particle 



exhibits intermittent behavior, which was first reported by Skotheim & Secomb (2007 1 . 



Figures [5] a-|6]c illustrate the complexity of the intermittent dynamics 

In order to compare with earlier works, we have considered an ellipsoid with excess 
area A = 0.02, which corresponds to a prolate ellipsoid with axes 0.9, 0.9 and 1. For 
this small excess area, the Keller-Skalak model (\ = 0) is in excellent agreement with 



the exact theory (5.2) and only slightly overestimates the viscosity ratio at the transition 



from TT to TB. However, when the elasticity is nonzero there is discrepancy between our 



results and those of Skotheim & Secomb ( 2007 ) , mainly because it is impossible to relate 



their phenomenological parameter U e to the elastic capillary number \- Our equations, 



however, agree with those derived by Kessler et al. (2009) (see Appendix Rl). 
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(b) 
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Figure 4. Effect of viscosity ratio on the swinging motion, (a) time evoiution of the inclination 
angle, (b) phase portrait. Parameters: A = 1 (red) , A = 10 (blue), and A = 20 (black). The 
excess area is A = 0.02 and elastic capillary number is \ — 10. 



5.1.1. Phase boundary 

In order to find the phase boundary between the SW and TB, the angle evolution 



equation (5.2 1 is rewritten as 



i> = -± + -^A(4>)ca8(2il> + 5(<f>)), 



where 



.4 



4A , 
1+ 30^ +4 * 



A 
30^ 



sin(20) 



1/2 



tan 



^Xcos(20) 



i+ySx^x 



(5.3) 



(5.4) 



At the saddle-node bifurcation from SW to TB, 

ipdt = V>(2tt) - ip(0) = . 



Since cos(2V> + S) is bounded between —1 and +1, the integrand, given by (5.3), would 
be zero if 

*A* = I 

VA 2 
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(a) 



(b) 



(c) 
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Figure 5. Intermittent dynamics: (a) time evolution of the inclination angle, (b) phase 
portrait, (c) zoom into the phase portrait. Parameters: A = 0.02, \ — 12.5, A = 20 . 



In order to estimate the critical value of the elastic capillary number \ c , we observe that 
the function A(x) has a minimum when 



sm(2(f>) = min { -2\j — X , -1 



(5.5) 
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(a) 



(b) 



(c) 
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Figure 6. Intermittent dynamics: (a) time evolution of the inclination angle, (b) phase 
portrait, (c) zoom into the phase portrait. Parameters: A = 0.02, \ — 21, A = 15 . 



At the bifurcation j4 mm = A* , which leads to (assuming that \ < *\/307r/4A) 



Xc 



30tt / A 1 
4AV ~J¥' 



(5.6) 



Note that the above derivation also indicates that this boundary should be an upper 
boundary. Figure [7] shows the phase boundary between SW and TB /intermittent dy- 
namics. 
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Figure 7. Critical elasticity number vs viscosity ratio. The points denote transition from SW 
to intermittent or T B behavior as determined from the return maps. The red line is the up- 
per boundary (5.61, Xc 1 — V 1 — A 2 , and the blue line is the Kessler et 
Xc 1 = 1 - A/2 ; the notation is A = v /A/30vr(23A + 32)/8 and x = 2xVA/30tt~ 



(2009 1 's result 



5.1.2. Weak elasticity 

In order to make further progress analytically, let us consider the limit of weak elas- 
ticity, x "C 1. In this case, the effect of shear elasticity is analyzed as a small per- 
turbation around the stationary fixed point corresponding to a tank-treading vesicle, 
tp(t) = ip(x = 0) + s(t), where ip(x = 0) = l/2arccos ( vA/2/il. After linearization of 
(5.2) we obtain a linear forced oscillator equation 



2ft 

ft 



2/j 



x- 



V307T 

It can be integrated to give for the swinging motion 



sin t . 



where 



/3 = 



s(t) = acos(i + (3) . 



16% 



y / (23A + 32) 2 



<30ir\ 2 



(5.7) 



(5. 



(5.9) 



Interestingly, the phase angle (3 does not depend on the elasticity. The amplitude of the 
oscillations increases linearly with the elasticity number, i.e., increases with increasing 
elasticity or decreasing shear rate. Also, at a given elasticity number Xi the oscillations 
decrease with viscosity ratio. 

Note that the above analysis is limited to viscosity ratios such that the inclination 
angle ip(x — 0) is not small, i.e. A ~ 0(1). Only in this case, during the linearization we 
can neglect a term of the type sx but retain the term ssin^(x = 0) . 

5.2. Effect of the deformability: analysis of the full system 

The deformability introduces another parameter to the phase digram - the initial condi- 
tion for the shape, ro, which measures the asphericity along the flow-vorticity direction. 
In the previous section, we saw that the initial condition g^o = and ro 7^ imply 
constant ro = \/A/2 and the shape is a fixed ellipsoid whose major axis lies in the plane 
perpendicular to the undisturbed vorticity. In the limit of high resistance to shearing 
X — > 00, the particle motion is given by the C = 00 Jeffrey orbit ( Jeffrey|[J923 1 . 
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2r /\/A 

Figure 8. Phase boundary between swinging and tumbling. Critical elasticity as a function of 
2ro/vA for different viscosity ratios A = (green), A = 5 (red diamonds), A = 10 (blue crosses), 
A = 20 (black circles). A = 0.02. The line is l/r 



If 5°o 7^ but r = the capsule attains a stationary ellipsoid shape. In the rigid 
ellipsoid limit, \ ~ * °°j the motion is equivalent to the C — Jeffery orbit (spinning 
around the axis of symmetry, which is parallel to the undisturbed vorticity ) ; the shape 
is axisymmetric and characterized by /20 = -\/A/2 and /2±2 = 0. For finite x the shape 
is a general ellipsoid. 

The most important consequence of the deformability is that it suppresses the inter- 
mittency. It appears that the intermittency is an artifact of the reduced model. Once the 



full set of equations (4.5) and (4.6 1 is solved, no intermittency is observed, only very long 



transients. This is in agreement with the numerical simulations by Kessler et al. (2008), 



and in contrast to the conclusion by Noguchi ( 2009 ) that deformability does not change 



qualitatively the dynamics. The latter analysis, however, ignores the effect of the tension 
due to the the area constraint. 

The phase diagram for a deformable particle is presented in Figures [8] and [9j Figure 
[8] shows that as r — > 0, the critical elasticity number diverges as l/r . Decreasing r , 
which allows for more deformation, widens the swinging region, as seen in Figure [9j 
Figure 10 illustrates the evolution of the limit cycle in the R — ip phase space. The 
decreasing radius with increasing tq indicates that enhanced deformability decreases 
the amplitude of the oscillations. Figure [Tl] shows that increasing the elasticity number 
(decreasing the shear rate) x increases the amplitude of the oscillations, and eventually 
causes tumbling. Finally, Figure [12] demonstrates that increasing the viscosity ratio A 
increases the amplitude of the oscillations. 



6. Conclusions 

We have considered the dynamics of a deformable membrane-encapsulated fluid par- 
ticle in steady shear flow. We have developed an analytical solution in the asymptotic 
case where the deformation is limited by small excess area. The theory accounts for 
the membrane area-incompressibility, resistance to shearing, and ellipsoidal unstressed 
shape. The solution of the evolution equation allowed us to construct the phase diagram 
and to derive analytical results for the phase boundaries. 

Our analysis clarifies the physical basis and limitations of earlier phenomenological 
models, which are restricted to fixed shape and compressible membrane. In particular, 
we show that the reported intermittency is an artifact of the shape-preservation. 
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Figure 9. Phase boundary between swinging and tumbling. Critical elasticity as a function of 
viscosity ratio for different deformability 2ro/\/A: 0.2 (green triangles) , 0.5 (red crosses), 0.8 
(blue diamonds), 1 (black circles). A = 0.02. 




2.R/VA 

Figure 10. Effect of r on the limit cycle at fixed \ = 20, A = 5, A = 0.02. 2r /V~A = 1 (red), 
0.8(blue), 0.5(green), 0.1 (black) and O.Ol(magenta) 




2.R/VA 

Figure 11. Effect of x on the limit cycle at given 2r /\/A = 0.5, A = 5, A = 0.02. \ : 
l(magenta), 10 (red), 30(blue), 50(green), 67(black). 
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-T0.0 0.2 0.4 0.6 0.8 1.0 3 

Figure 12. Effect of A on the limit cycle at given 2r /\/A = 0.5, \ = 20, A = 0.02. A 
0(magenta),10 (red), 20 (blue), 30(green), 40(orange), 48(black). 
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Appendix A. Relation to the analysis of Kessler et al. (2009) 



In Kessler et al.'s paper Kessler et al. ( |2009 ), angles E and A (not to be confused by 
the excess area in our notation) were employed. They are related to the inclination angle 
ip and the tank-trading angle cf> as follows 



E = -0-J, A = <£-2V + J. 
Starting from their equations (21) and (22) 
<9 r E=A 
d T A =A + 4sinEsinA + 2(x~ 1 - l) cos (E - A) 



(Al) 



(A 2) 



we recover our equations (4.5) and (4.6), if we identify A and \ 1 from their notation 
with 



\& 



A 23A + 32 „_! 
307T 8^ ,X X V 30,.- 



A 



(A3) 



Performing the same analysis on (A2) as in section 5.1.1 yields the phase boundary 
reported by Kessler et al., namely 



Xc 



V30vr 



1 ^-^(23A + 32) 

16730^ 



(A 4) 



Appendix B. Shape evolution equations for the deformation of an 
ellipsoidal particle in shear flow 

Ellipsoidal deformation is characterized by only j ' = 2 modes. Let us introduce fj m 

fjm + ifjrn and fj-m = (-l) m (/jm - tf'j m )'> Same for tllC 9jm 
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The shape evolution is described by 
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a/20 A h » 4h 
^ r =4^/ 20 / 22 +X 



9/2 



AVM/f 

4/i 



ff 2o (A - 2/f ) - 4/ 20 (322/22 + .922/22) ] (B 1) 



dt /22+4 A /22/22+X Ax^0tt 



322 ( A ~ 4 (/22J J - 2/22 (.920/20 + 23 22 / 22 

(B2) 



Of: 



dt - h+ ^ +i A [H +X AVMn 



Ah 



922 A-4 /. 



at 

Next we describe the solution. 



<9fl 22 " dg 22 > dg 20 n 

" " W 322 -^T = w #22 ~^t- = 



(•;/ 



9* 



2/22 (320/20 + 2g 22 / 22 
(B3) 

(B4) 



Appendix C. Formalism 

C.l. Velocity fields and hydrodynamics stresses 
Veloci ty fields are described using basis sets of fundamental solutions of the Stoke s equa- 



tions (Schmitz & Felderhof 
defined in Appendix [F] 



1982 



Cichocki et al. 



1988 



Vlahovska et al. 



2009a I, u 



jrnq' 



V° X (r) = J2 C Tmg [ U tmq( r ) ~ U jraq( r )] + ^ C i^ U «( F ) > 

v in (r) = ^ c J m 9 u J + ™?( r )- 



(Clo) 
(C16) 



On a sphere u+ (r = 1) = u7 (r = 1) and the velocity fields given by (CI I arc 
continuous. 

The hydrodynamic tractions exerted on a surface with a normal vector n are n • T. 



They are expanded in vector spherical harmonics ( E 3 ) 



t — n • ± ^jrnqYjr 



(C2) 



In the particular case of a sphere characterized with a normal vector f , the amplitudes 
of the viscous tractions and the velocity field are linearly related 



tjmq = Z^ c J m i' \ q'q ~ q'q) + 2_^ c J m <i'®q 



(C3o) 



t in -Vc- /9+, 
q' 



(C36) 



where 8i are obtained from the velocity fields (F1|-(F2| ( Blawzdziewicz ei a2.|2000), 

(C4) 
2j + l + f y 



©^ (i) 



2.7 + 1 




-3(t 

i - 1 0' 

^(^ n > i.-< 
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e;,' V) 



I -2j-l 

-j - 2 



\<m) 



o 



3 bii 



-2.7 - 1 



1 \ 



(C5) 



7 + 1 / 



C.2. Boundary conditions: Stress balance 
The stress balance in terms of spherical harmonics reads 



jmq ~ Z jmq ' 



(C6) 



Tangential stresses correspond to the q — 0, 1 components, and the normal stresses - to 



"jmq 



( - /Cl k, V'jmq • ''jmq) ' X^jmq ' 



q = 2. The hydrodynamic tractions are given by (C3&I. The membrane tractions are 
flVlahovska fc Gracia|[2007l |Seifert|[l999| 

(C7) 

(C8) 

/; = -v'./(/ + lV)-;,„. (C9) 



where the bending stresses are 

t- m2 = JU + 1) U - 1) (i + 2) /,-„, , i? m0 

and the stresses due to membrane tension are 



t%n2 = 2(T jm + (To (j - 1) (.? + 2) f jr 



jmO 



Vj'O' + iV. 



jm • 



Note that the membrane tension is non-uniform under non-equilibrium conditions 



0~ — 0~Q + 2_^ 0~j m Yjm ■ 
jm 



The surface elastic stresses have only in-plane shearing component 

tf m , = , t» =2( 3 - 1) (j + 2) [j(j + l)]- 1/2 (f 3m - 9j? 



jm2 



(CIO) 



(Cll) 



where g 3m denotes the reference configuration. Membrane stresses do not involve a q = 1 
component. 



Appendix D. Solution 

Simple shear flow is defined as 
which translates into 



D.l. External flow 



v°°=yx 



2 j 2 

V°° = V V V C°° U+ 

Z_^i Z_^ Z_^, 3-mq jmq 

j=l m=—j 9=0 



2±20 — "FU/ Z i C 2 ±22 — TM 



i °101 



(Dl) 
(D2) 
(D3) 



D.2. Rotation: tank-treading frequency 
First, we find the velocity field amplitude Cj m \ using the tangential stress balance with 
q=l 

tf ml -\tf ml =Q. (D4) 



This gives 



Swinging of red blood cells in simple shear 
2j + l 



Cjml — c j m i 



2 + j + X(j-l) 
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(D5) 



The imposed shear flow ( D 3 ) has a rotational component only with j = 1 and hence the 
above relation reduces to 

doi=c? 01 (D6) 

i.e., the rigid body rotation is unperturbed by the particle, and the particle rotates with 
the same rate as the flow. 

D.3. Membrane incompressibility 

The local area conservation implies that the velocity field at the interface is solenoidal 
flSeifert|[l999l ) 

V s • v = at r = 1 (D 7) 

Therefore the amplitudes of the velocity field (CI) are related 

2 



GjmO 



vJU+V 

D.4. Tension 



Cjm2 



(D» 



The non- uniform part of the membrane tension, <7j m , is determined from the tangential 
component of the stress balance (C6), q = 0, 



°~jm — ka K 



oo 2(l+2j) , ^,00 3(2j+l) 



(2+j + AQ-l)) (j-lXj+2) 



(D9) 
D.5. Normal velocity and shape evolution 



The tension ( D 9 ) is substituted into the normal component of the stress balance ( C 6 ) , 
q = 2, to obtain the normal velocity c jm 2 

c,- m2 = C iro -(j+2)(i-l)d(A, j) _1 [<7a -1 C7(3 + l)(j(j + 1) + ffo))/jm + Wjm - ffim)] , 

(D10) 
where 



C im = d(A, j)- 1 cf mQ y/j(J+l) (2j + 1) + c- 2 (4j 3 + 6j 2 - 4j - 3) (D 11) 



and 



d(X,j) = (4 + 3j 2 + 2j 3 ) + (-5 + 'if + 2f)\ . (D 12) 

Finally, the motion of the interface is determined from the kinematic condition 



c)t 



Cj,n2 + v s ■ V/ at r = 1 



where v s = v cx (r = 1) = v ln (r = 1). At leading order, it takes the form 

dfjr 



01 



= c jm2 +iuj 7 -§f jm atr = l. 



(D13) 



(D14) 



U) = 1 is the local rate of rotation. Note that jL — iwm/2 is in fact the Jaumann derivative. 
Substituting c,- m 2 hi (D 13) yields the evolution equation for the shape parameters 



9f. 



^ = iu^f jm + C jm + Ca- 1 (Ti + r7 r 2 ) f jm + xT 3 (f jm - 9 jm ) + O (A) , (D 15) 
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where 

Ti = -U + m 1) (UU + l)] 2 ) d(X, J)' 1 , (D 16) 

T 2 = -(j + 2)(j - l)j(j + l)d(XJ)- 1 , (D 17) 

T 3 = -4(j + 2)(j-l)d(\,j)-\ (D18) 

The reference configuration simply rotates with the flow 



dg 



3 ,n — ;, ,m. 



<9f 



(D19) 



D.6. Isotropic tension 



The normal velocity (D10) and the shape evolution (D15l include the yet unknown 



isotropic membrane tension. It is expressed in terms of the shape modes and other known 
parameters in the problem using the area constraint ( Vlahovska fc Gracia| [2007) 

(i + 2)(j-i) 



A = Yl a U)fjmfj m , a(j) = 



(D20) 



Since A = 0, it follows that 



j:-(J)(f;J-t + fJ-§ 



= 0. 



(D21) 



Using (D 15 ) leads to 
Ca K 



o-o 



Yjjm a U) T 2fj m fj m 



^2o,(j) [C Jm f* m + Ca K ^i/j™/^ + X^3(fjm - gjrn)fj m ] ■ 

(D22) 



jm 



In order to clarify the physical significance of the isotropic tension, let us consider the 



particular case when only the ellipsoidal deformation modes, j = 2, are present. (D9| 
simplifies to 



.Ca K /lOvr u 2 

o-q = -6 - i— d — (/ 22 - /2-2J + ^xCa K 



We see that the tension varies with deformation. 



T + -r- (320/20 + 322/2-2 + .92-2/22) 



(D23) 



modes yields (B1)-(B4) 



Inserting the tension (D 22 ) into the shape evolution (D 15 1 and keeping only the j = 2 



Appendix E. Spherical harmonics 

A good reference on spherical harmonics is Varshalovich et al. (1988). The normalized 
spherical scalar harmonics are defined as 



Y jm (9, ip) 



2j + l (j-m)\ 
4-7T (j+m)! 



-l) m P™(cos(9)e imv , 



(El) 



where f = r/r, (r,9,(p) are the spherical coordinates, and P m (cos#) are the Legendre 
polynomials. For example 



^io = cos 9 . 

V47T 



(E2) 
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The vector spherical harmonics are defined as Blawzdziewicz et al. ( 2000 ) 

_ i 
yjmo = [i U + 1)] 2 rVaY jm , 

yjm2 = rYjm , ( E 3 ) 

Yjml = -if X VjmO 

where Vq denotes the angular part of the gradient operator. 



Appendix F. Fundamental set of velocity fields 

We reproduc 
list the express 
is described by 



We reproduce the velocity and stress fields from Vlahovska & Gracia ( 2007 ) . First we 

± 

jmq 



list the expressions for the functions uf m „ (r, 0,p). The velocity field outside the vesicle 



u" m0 = \r~i (2 - j + jr- 2 ) y jm0 + \r~i [j (j + 1)] 1/J (l - r~ 2 ) y jm2 , (F la) 

1^1=^-^;™!, (F16) 



jml 

1/2 



U jm2 ~ 2' 



\r- ] (2 - j) (i^) (1 - r- 2 ) y jm0 + \r^ {j + (2 - j>- 2 ) y jm2 . 

(Flc) 
The velocity field inside the vesicle is described by 



u 



r J - X (~(j + 1) + (j + 3)r 2 ) y jm0 - \r^ [j (j + 1)] 1/2 (l - r 2 ) y jm2 , 

(F2o) 
< rol =Wy im i, (F26) 

W>- o- (3+j)(^) 1/2 (l-r 2 ) yjm0 + i^- 1 ( J +3-( J + l)r 2 ) yjm2 . 

(F2c) 
On a sphere r = 1 these velocity fields reduce to the vector spherical harmonics defined 



jmO 2 



u+ = Ir?- 1 



by (E3I 



U jm ? — yjmq ■ (F 3) 

Hence, v.^ m0 and u^ ml are tangential, and u - 2 is normal to a sphere. In addition, u ■ 
defines an irrotational velocity field. 
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